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We construct an exactly soluble Hamiltonian on the D=3 cubic lattice, whose ground state is a 
topological phase of bosons protected by time reversal symmetry, i.e a symmetry protected topo- 
logical (SPT) phase. In this model anyonic excitations are shown to exist at the surface but not in 
the bulk. The statistics of these surface anyons is explicitly computed and shown to be identical 
to the 3-fermion Z2 model, a variant of Z2 topological order which cannot be realized in a purely 
D=2 system with time reversal symmetry. Thus the model realizes a novel surface termination for 
SPT phases, which only becomes available in D=3: that of a fully symmetric gapped surface with 
topological order. The 3D phase found here is also outside the group cohomology classification 
that appears to capture all SPT phases in lower dimensions. It is identified with a phase previ- 
ously predicted from a field theoretic analysis, whose surface is roughly 'half of a Kitaev Eg phase. 
Our construction utilizes the Walker- Wang prescription to create a 3D confined phase with surface 
anyons, which can be extended to other symmetry classes and topological orders. 



Introduction Recently, there has been much inter- 
est in studying gapped phases of interacting bosons that, 
despite having only conventional gapped excitations in 
the bulk, are topologically distinct from the trivial phase. 
Such short range entangled (SRE) phases are signifi- 
cantly simpler than topologically ordered phases, such 
as fractional quantum Hall states and gapped spin liq- 
uids, whose bulk anyonic excitations reflect their long 
range entangled naturepQ [2]. Nevertheless, SRE topo- 
logical phases display a host of interesting phenomena: 
for example, like the noninteracting fermionic topolog- 
ical insulators, they possess exotic surface states [3H5]. 
However, in contrast to the fermionic case, topological 
phases of bosons cannot be realized without interactions. 
Often, SRE topological phases are well defined only in 
the presence of a symmetry. In those cases, they are 
termed Symmetry Protected Topological (SPT) phases. 
The best known example is the Haldane antiferromag- 
netic chain of S=l moments in 1+1D, which is in a SPT 
phase protected by spin rotation symmetry [7]. Bosonic 
topological phases in 1+1 dimension occur only in the 
presence of a symmetry group G, and are rigorously ar- 
gued to be classified by the second cohomology group 

Recently it has been realized that analogous states ex- 
ist in higher dimensions [T3HT7] , and could potentially be 
realized as ground states of frustrated magnetic insula- 
tors or ultra cold systems of lattice bosons [T8H2T] . In 
2+ ID the situation is especially rich, since SRE topo- 
logical phases that exist independent of symmetry are 
possible. Bosonic phases with c_ = 8n chiral bosonic 
edge modes (with n any integer) were shown by Ki- 
taev to be compatible with SRE. This leads to a series 
of topological phases with quantized thermal Hall con- 
ductivity K xy jT = 8n[22 [23 . Thus different integers n 



lead to physically distinct transport properties, giving 
a sharp criterion that can be used to distinguish these 
phases. We will refer to the n = 1 member of this se- 
quence as the Kitaev Eg state. In the presence of sym- 
metry, more topological phases have been predicted, and 
it was proposed that SPT phases are classified by the 
higher dimensional cohomology groups H d+1 (G, U(l)) in 
d + 1 dimensions fl3\ [14]. The predicted set of phases 
from this group cohomology approach was reproduced 
using other techniques in a number of specific cases: 
in 2D via dualities [15] and the Chern Simons K-matrix 
theory[16, 24 , and by various field theoretical approaches 
in 3D [HI 125]. 

Intriguingly, however, the field theoretical 
approach [17 predicted new SPT phases in 3+1 D 
protected by time reversal symmetry T. Thus, in 
contrast to the single nontrivial phase predicted by 
group cohomology, two distinct types of nontrivial 
phases were obtained. One of these nontrivial phases 
has a close connection to the 2D Kitaev Eg state: a 
domain wall between opposite T breaking regions on 
the surface carries a chiral central charge c_ = 8, i.e. 8 
chiral boson modes, and corresponds to the edge of the 
Kitaev Eg state. Thus each domain may be associated 
with half this number of chiral edge states, a feature 
that is forbidden in a purely 2D system in the absence 
of topological order. This bosonic state is analogous to 
the 3+1D free fermion topological superconductor (class 
Dili), where a surface domain wall between regions of 
opposite T breaking is associated with a single chiral 
Majorana mode. Hence we will refer to the bosonic 
phase as the 3D Bosonic Topological Superconductor 
(BTSc). This state was proposed to be outside the 
group cohomology classification, and although it was 
discussed as a physical possibility in Ref. [17], more 
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direct evidence for its existence would be welcome. 
This would shed light on the exact relation between 
group cohomology and SPT phases. Here we irrefutably 
demonstrate its existence by providing an exactly soluble 
model that realizes this phase. 

Surface states of 2+ ID SPT phases must either be gap- 
less or spontaneously break symmetry. This is analogous 
to the surface states of free fermion topological insulators 
which are required either to be metallic or break one of 
the protecting symmetries. On the 2D surface of a 3D 
SPT phase, a novel possibility may arise - a fully gapped 
and symmetric state is allowed if the surface develops 
topological order [T71I26]. However, the resulting topolog- 
ically ordered surface is anomalous - i.e. it implements 
the global symmetries in a way that cannot be realized 
in a strictly 2+1D phase [T7]. In the context of the 3D 
BTSc, the topologically ordered surface is predicted [T7] 
to realize the "3-fermion Z2 state" , with T. Recall that 
the regular Z 2 topological order, as embodied by the toric 
code or a Z2 gauge theory, has three nontrivial anyons: 
the electric (e) and magnetic (m) charges and their bound 
state (e). e and m are bosons, and e is a fermion; all ex- 
citations have mutual tt statistics. In the 3-fermion Z2 
state, all three nontrivial anyons are fermions, again with 
mutual 7r statistics. ReflTTl proposed that this is a possi- 
ble surface termination of the 3D BTSc with time rever- 
sal symmetry. A strictly 2D realization of the 3-fermion 
Z2 state always breaks time reversal symmetry since it 
is associated with c_ = 4 chiral edge boson modes [22] . 
Since this surface state in the presence of T can never 
be realized as a 2D system, it qualifies as the protected 
surface state of a 3+1D SPT phase. Our exact solution 
produces precisely the 3-fermion Z2 state in the presence 
of a boundary, while preserving T. Thus not only does 
our solvable model exhibit an interesting bulk phase; it 
also turns out to have the most exotic possible surface 
state. The key observation is that although a 2D real- 
ization of the 3-fermion Z2 phase always breaks T, the 
self and mutual statistics of the anyons only involve real 
±1 phases, allowing us to construct a T symmetric 3D 
model. 

We use the Walker- Wang [23 [28] prescription for con- 
structing an exactly soluble model of a 3+ ID phase 
given a 2+1 dimensional topological quantum field the- 
ory (TQFT)[29]. Let us assume that the anyons of the 
topological theory are such that the only particle that 
has trivial mutual statistics with all the others is the 
identity particle (a Modular Tensor Category). In this 
case the Walker- Wang models realize a confined 3+ ID 
phase [27j [28]. However, the surfaces are naturally en- 
dowed with the topological order of the 2+1D TQFT. 
Since the bulk is not topologically ordered, and, in the 
absence of symmetry, there is no obstruction to destroy- 
ing the topological order at the surface, these models lead 
to trivial insulators. [30] This is an ideal starting point for 
constructing SPT phases, which, in the absence of sym- 



metry, are connected to the trivial insulator. As we will 
see, in the presence of a global symmetry, the Walker- 
Wang models can have topologically protected surface 
states. 

In the remainder of this paper, we first discuss some 
background about the 3-fermion Z2 model and the 
Walker- Wang construction. We explain why in these 
models the bulk is confined, but deconfined excitations 
of a specific kind exist at the boundary. Since these ar- 
guments are rather general and intuitive we review them 
before presenting a microscopic model. We then discuss a 
3D model on the cubic lattice based on the 3-fermion Z2 
state, which is composed of a sum of commuting projec- 
tors and hence exactly soluble. The model is shown to be 
T symmetric, confined in the bulk but with deconfined 
excitations at the surface, which are explicitly shown to 
have the same statistics as the 3-fermion Z2 state. 

Background and Walker- Wang Construction 

We begin by describing the 3 fermion state in 2D, intro- 
duced by Ref. [22]. We may think of the three types of 
fermions as fermionic Z2 charges (e), fermionic Z2 fluxes 
(m), and a bound state of a charge and a flux (e). Be- 
cause the charges acquire a tt Berry phase upon encircling 
the Z2 fluxes, it can be checked that e is also a fermion. 
Moreover, the three species of anyonic excitations all have 
mutual semionic statistics, i.e. braiding one around an- 
other induces a phase factor of —1. Due to the symmetry 
in these statistics, we call this topological state the 'three 
fermion Z2 model'; in practice it is irrelevant which one 
of the labels {e, m, e} we assign to the flux and which to 
the 'original' fermionic charge. We will use the label 1 to 
designate the vacuum. 

One important property of this three fermion model 
in 2D is that it breaks time reversal symmetry. This be- 
comes evident when the system is put on a disc whose 
boundary will have chiral edge modes. In particular, us- 
ing the relation between the topological properties of the 
bulk anyons and the edge chiral central charge 



a V a 



it is easy to check that c_ = 4 (mod 8) , with quantum 
dimensions d a = 1, and topological spin a = — 1 for e,m 
and £, corresponding to fermions. Therefore, the bulk 
topological property dictates that this state can not be 
realized in a 2D time reversal invariant system. 

An explicit field theory of this 2D state can be written 
using an Abelian Chern Simons theory with four U(l) 
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gauge fields: 
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This is the Cartan matrix of SO (8). Note that the inverse 
matrix is: 
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which clearly demonstrates the mutual statistics of the 3 
fermion model, obtained by inner products Oij = 2nlf • 
K~ x • lj , while self statistics is given by: 6i = nlf • K~ x • k 
where k are integer vectors representing the quasiparti- 
cles. The eigenvalues of K so ^ are all positive, indicat- 
ing that all four edge modes are chiral (propagate in the 
same direction). Therefore the 3 fermion state explicitly 
breaks time reversal symmetry, when realized in 2D. 

To realize such a surface state, we will use a Hamil- 
tonian based on the Walker- Wang [27] [28] construction. 
This leads to a ground state wave function which is a su- 
perposition of loops (more precisely "string nets" , in the 
sense of [29]) labelled by the anyon types. The amplitude 
for a given configuration of these loops C in the 3+ ID 
wave function ^3b(C) is determined by the expectation 
value of the corresponding Wilson loop operators in the 
2+1D TQFT; i.e. we set: 

^3D(C) = (Jy(C)) 2 +lTQFT 

This is similar in spirit to eg. Quantum Hall wavefunc- 
tions, which are related to the space-time correlations of 
their edge states. Here, since we demand a topologically 
ordered boundary state, the expectation values are taken 
in the boundary TQFT. Below we will argue more phys- 
ically why the bulk wave function encodes the statistical 
interactions of the surface anyons, but why the bulk is 
itself not topologically ordered. 

The wave function hence contains three different colors 
of loops corresponding to the three species of fermions. 
Any two colors can merge into the third, in accordance 
with the fusion rules of the TQFT. Each loop configu- 
ration comes with a specific phase due to the twisting 
and intertwining of the fermion world-lines. Since this 
twisting may depend on the angle of view, in order to 
calculate the phase, we need to pick a particular projec- 
tion of the 3D wolrd-lines onto 2D ones. The projection 
we will use is shown in Figjl] Having fixed a projection, 



the phase factor can be obtained using the braiding rules 
(as shown in Fig{T| given by the R matrix: [31] 
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Isolated loops can shrink to the vacuum without an ex- 
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FIG. 1. (color online) Braiding rules for strings in the ground 
state wave funciton. (a) applies to strings of the same color 
while (b) applies to strings of different colors. 

tra phase factor (the quantum dimension of our fermions 
is 1). The ground state wave function is a superposition 
of all allowed loop configurations weighted by the corre- 
sponding phase factor. When the system has a surface, 
the same graphical rules can be used to determine the 
wave function for the loop gas after fixing the projection. 

An important feature of this state is that, because the 
braiding rules involve no complex numbers, the ground 
state wave function is real, hence symmetric under the 
time reversal operator T that acts by complex conju- 
gation. This is even true for the wave function on a 
3D manifold with boundary. Therefore, no time reversal 
symmetry breaking occurs either in the bulk or on the 
surface. 

On the other hand, the graphical rules allow us to see 
that the surface of the 3D model has the anyon content 
of a state which would violate time reversal in 2D. In 
particular, anyonic excitations can be created by adding 
open strings to the surface. The wave function becomes a 
superposition of all string configurations in which strings 
end at the positions of the excitations. Then we can 
check the statistics of the excitations by tracking these 
open strings. Suppose we exchange two string ends of 
the same color (say red, as show in Figj2](a)) by crossing 
two red string segments to the surface. (Fig{2|a) shows 
one possible configuration.) This twist in the string con- 
figuration (relative to the string configuration before ex- 
change) can be removed to bring the strings back to their 
original form, but this results in a (— ) sign. Therefore, 
exchanging end of strings of the same color adds a (— ) 
sign to the total wave function, which is equivalent to 
saying that the ends of the strings are fermions. Sim- 
ilarly one can check that string ends of different colors 
have mutual semionic statistics by braiding them with 
linked loops on the surface as shown in Figj2] 

Open strings in the bulk also create excitations in the 
ground state. However, the excitation energy grows lin- 
early with the string length, leading to confinement of the 
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(a) 



(b) 



FIG. 2. (color online) The anyonic excitations on the surface 
are created by open strings. At the ends of the strings are 
three species of fermions (corresponding to three colors of 
open strings) which have mutual semionic statistics. This can 
be seen from the braiding statistics of the strings generating 
(a) the exchange and (b) the braiding of the end of strings. 



particles at the ends of the strings. To see the confine- 
ment, consider an open string (for example blue) in the 
bulk which is circled by a small ring of a different color 
(for example red), as shown in Fig J3] Using the braiding 





FIG. 3. (color online) Open strings in the bulk change the 
quantum fluctuation phase factors of loops along its length, 
which costs finite energy. Therefore, the end of strings are 
confined in the bulk. 

rule between the loops we find that the linking between 
the ring and the open string can be removed together 
with a (— ) sign. That is, the open string changes the 
quantum fluctuation phase factors of small loops along its 
length, which costs finite energy. Therefore, the string's 
endpoints cannot be separated very far, and the fermionic 
excitations in the bulk are confined. In fact, any set of 
strings with consistent braiding and fusion rules given by 
a unitary braided fusion category can be used to write 3D 
string-net wave functions in a similar fashion. As shown 
in Ref.[28 , the bulk of the state has no deconfmed ex- 
citations, and hence no nontrivial topological order, as 
long as each string has nontrivial statistics with at least 
one of the strings. The corresponding category is said to 
be 'modular'. However, open strings lying on the surface, 
where no loops can encircle them, may give rise to decon- 
fined excitations. As we prove in more detail later, for 
the particular case of the 3-fermion Z2 model the excita- 
tions at the ends of open strings are indeed deconfmed 
on the surface, giving rise to three species of deconfmed 
fermions with mutual semionic statistics on the 2D sur- 
face of a 3D confined model, as promised. 

The surface of the 3D model hence realizes a state 
which is impossible in pure 2D systems with time re- 
versal symmetry. Such a surface state is an indication 
of the nontrivial SPT order in the bulk and cannot be 



removed wihtout breaking time reversal symmetry or go- 
ing through a bulk phase transition. To see this, we 
can imagine destroying the topological order by adding 
a 2D surface layer with the same types of excitations 
and binding the corresponding pairs of excitations to- 
gether; however, as this 2D system has chiral fermionic 
edge modes [22 , this operation must break time rever- 
sal symmetry. In fact, if we perform the above surface 
confinement operation using opposite T breaking states 
in different halves of the sample, this creates a domain 
wall. Since each 2D state has edge chiral central charge 
c_ = 4, the combination leads to a net chiral central 
charge c_ = 8 along the domain wall. This is expected 
for the non fractionalized surface of the 3D Bosonic Topo- 
logical Superconductor discussed in [17] . It also indicates 
that when we destroy the topological order of the surface, 
this phase breaks T symmetry and displays a thermal 
analog of the quantized magneto-electric effect [3H5]. 

One may directly argue that realizing a topologically 
ordered phase which transforms under symmetry in a way 
this is forbidden in 2D immediately leads to a protected 
surface state. Assume for the moment that the opposite 
is true - i.e. there is some way to confine the 3-fermion 
Z2 state without breaking T symmetry. Then, one can 
make a slab of the 3D phase with well separated top and 
bottom surfaces, and eliminate the surface state on the 
bottom side. Now, consider shrinking the slab till the 
2D limit is reached. Since the bottom side and bulk are 
gapped, it should be possible to retain the original surface 
state on the top surface, without changing the symme- 
try. Now we have achieved a contradiction, by having 
produced a 2D realization of an 'impossible' 2D state. 
Therefore our assumption that it is possible to eliminate 
the surface state without breaking the symmetry must 
have been incorrect. 

Exactly Soluble Model of a 3D SPT Phase with 
Surface Topological Order We now describe an ex- 
actly solvable 3D model which realizes these properties. 
The Hamiltonian we consider is defined on the cubic lat- 
tice. The Hilbert space consists of 4 possible states on 
each edge, which we call |1), |e), |ra), and \e). The Hamil- 
tonian is a sum of commuting vertex (V) and plaquette 
(P) terms: 



(6) 



Here the first sum is over all vertices V of the cubic lat- 
tice, and Ay gives an energy penalty to certain combi- 
nations of the possible states of the 6 edges adjoining 
V. The combinations that are energetically cheap are 
dictated by the "fusion rules" for the 3-fermion model: 
in the ground state, combining the anyons on all edges 
adjoining V must give the identity. Thus closed loops 
of a single anyon type are allowed, as are vertices where 
the three nontrivial anyons fuse together. It is useful to 
think of the states |e) and \m) as carrying the non-trivial 
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charge of two Z2 symmetries (Z^ and Z^), and the \e) 
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state as charged under both Z^ and Z^, while |1) car- 
ries no charge. The vertex term Ay then simply enforces 
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the conservation of Z^ and Zi>' " ; charges at V. Some of 
the configurations that satisfy this condition are shown 
in Fig. [4] 

The plaquette term is a little more involved. First, 
define operators 5^ by 

5"|1) = |m), 5» = VI1> + I^a||A) (7) 

where (/i, z/, A G {e, m, £}), with |e Miy A | = 1 when /i, z/, and 
A are some permutation of e, m, and e, and otherwise. 
We can think of as raising the appropriate Z2 charge 
carried by the edge upon which it acts. For example, 
S e adds a non-trivial Z^ charge to the link, changing 
an \e) edge into an edge that carries trivial charge under 
both symmetries (i.e. the |1) edge), the \m) edge into 
an edge that is charged under both symmetries (the \e) 
edge), and so on. 

The plaquette operator Bp acts on <9P, the set of edges 
bordering P, and two other special edges O and U. To 
define these, we fix a 2d projection of our 3d lattice once 
and for all, and let O and U be the two edges which 
adjoin a vertex of the plaquette P and project inside it 
under the 2d projection (Fig. [5|. O is the edge which 
sits above (or "over") the plaquette, and U the one that 
sits under it. Let \a) and \ j3) be the states sitting on the 
O and U edges respectively. [32 Then 
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Here 6° and G u are 4 by 4 matrices of signs, indexed by 
the labels 1, e, m, e, and a, f3 are indices corresponding to 
the states on the O and U edges respectively. These are 
related to the R matrices introduced in Eq. ([5| by 
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FIG. 4. (color online) The Low energy Hilbert space of 
the lattice model consists of loops of three colors that satisfy 
fusion rules at the vertices - i.e. they are either closed loops 
of a single color, or segments of three colors can meet at a 
vertex. 
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FIG. 5. Choice of edges on which Op acts, for the three 
different types of plaquettes in the lattice. In the chosen pro- 
jection O edges (red) cross over the plaquette P, and U edges 
(blue) cross under it. 



Explicitly, the e, m, e parts are given by 
9^ = 




(10) 



and 0° = (O* 7 ) . Since the matrix elements of H are 
real, the Hamiltonian is invariant under time reversal T, 
where T is defined to be complex conjugation of the many 
body wavefunction in our l,e, m, e basis. Note that this 
time reversal operator satisfies T 2 = 1. 

With the definition (J8|, all of the terms in the Hamil- 
tonian commute. Since Bp raises the Z2 charge on a pair 
of edges at each vertex, it conserves both Z2 charges, and 
thus commutes with Ay, at every vertex V. To see that 
[Pp 1 ,Pp 2 ] = 0, we must check that the signs of all ma- 
trix elements are independent of the order in which the 
plaquette terms act. If Pi contains no O or U edges of 
P2 and vice versa, this follows immediately. Further, if 
Bp x raises the U edge of P2, then Bp 2 also raises the O 
edge of Pi (and vice versa). In this case the phase terms 
which could potentially depend on the operator ordering 
are: 
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where [12 x i indicates the state on edge i after acting with 
S^ 2 . One can check using Eq. (10) that these are equal 



for every choice of /ii, /i2? and so that [Bp 1 ,Bp 2 ] = 
0. 

Because all operators in the Hamiltonian commute, we 
can solve for the spectrum exactly. Before turning to this, 
let us develop further the algebraic formulation suggested 
by the mapping to Ising charges Z^ 6 '^ . Let the Pauli ma- 
trix af (rf ) measure the Ising charge Z^ (Zr, 77 ^) flowing 
through the link i. Then, the vertex term can be written 
as: 
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where the product is taken over the six links emanating 
from a vertex. The plaquette operators in Eq. ^ 
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are a product of three operators S^, 0^ and 0j^, where 
the latter two are interpreted as operators diagonal in 
the a x , r x basis and given by the matrix elements in 
Equation [lO] These operators are nontrivial only for \i = 
{e, m, e}, so we display them as a 3 component column 
vectors: 
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It is readily seen that the different plaquette operators 
Bp commute with one another: for any two different 
links 1, 2, 0^2 commutes with S%Q®i for all pairs 
/i, v. 

Excitations in the Bulk: We now show that there are 
no low energy point-like excitations in the bulk, by es- 
tablishing all bulk excitations are confined. This points 
to a short range entangled bulk state without topological 
entanglement entropy or 3 dimensional topological order. 
Essentially, the string operator that creates point excita- 
tions also creates a string of excited plaquettes, leading 
to linear confinement of the ends of the string [28 . We 
will also show that deconfined excitations exist on the 
surface. 

To create a pair of point particles we consider an oper- 
ator that violates the vertex term Ay at vertices V\ and 
V2. The simplest possibility is to act with S- 1 on all edges 
in some path C12 that connects V\ to V2 (see Fig. [6|. For 
example, we can attempt to create a pair of e particles 
by acting with n^ec^ a i • However, in doing so we nec- 
essarily change the state on the O and U edges of certain 
plaquettes along this path, creating a pair of strings of 
defective plaquettes running parallel to C12. For a path 
in the y direction, affected plaquettes are those in the xz 
plane that lie above and to the right (U edges) or be- 
low and to the left (O edges) of the path. Therefore the 
energy cost to separate the ends of the string operator 
grows linearly with distance. We can reduce the number 
of affected plaquettes per unit separation from 2 to 1 by 
using the modified string operator: 
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where *Ci2 is the set of edges that are crossed by a curve 
that runs parallel to C12, but is displaced slightly in the 
—x + y + z direction, and *(^^ er / under are subsets of 
edges that cross over (under) the path C12 in our projec- 
tion (colored purple and green respectively in Figure [6| . 



The terms in S M exactly cancel the phases that lead 
to non-commutativity with the plaquettes whose O and 
U edges are raised; however, the new operator fails to 
commute with the string of plaquettes (shaded blue pla- 
quettes in Figj6| threaded by the curve used to determine 
*Ci2 (dashed blue line in Figj6|. For example, to create 
a pair of e particles at Vi^ we act with S e = a z along 
the bonds connecting Vi ? 2, but also with the operators 
0^ = a x and 0^ = a x r x on the purple and green bonds 
respectively in Figure [6j It can be explicitly checked that 
this string operator does not commute with the shaded 
plaquettes but does commute with all other plaquettes. 
It is not possible to further reduce the number of violated 
plaquettes for a given Ci2[27l[28]. Observe that if we ter- 
minate the system just above the plane of the curve C12, 
the defective plaquettes would not be included in the 3D 
lattice, leading to deconfined surface excitations. This is 
indeed the case as we will see below. In fact, the sur- 
face has three species of deconfined fermionic excitations 
with mutually semionic statistics, corresponding exactly 
to the excitations in the 3 fermion model. 




FIG. 6. Excitations in the bulk are confined. The path C12 
is shown in red; the displaced path used to determine *Ci2 
is indicated with a dashed blue line. Edges that cross un- 
der (over) this path are colored green (purple). The violated 
plaquettes (shaded blue) are those that are threaded by the 
dashed blue line. 

Deconfined Surface Excitations: Consider terminating 
the 3D system with horizontal plaquettes in the x-y plane 
as shown in Figure [7[ The Hamiltonian for the surface 
vertices and plaquettes is derived from the bulk Hamil- 
tonian by simply retaining terms for the existing bonds 
on the surface. Thus both the horizontal plaquette and 
vertex terms at the surface each involve a product of 
five bond operators. With these boundary conditions the 
Hamiltonian is again a sum of commuting terms, and we 
can determine its spectrum exactly. 

Deconfinement occurs at the surface because, if C12 
lies entirely on the boundary, acting with an operator of 



the form (15) creates a pair of vertex defects but com- 



mutes with all plaquette projectors. (Note we may omit 
the 0° operators, since those bonds are absent). Thus, 
creating a pair of e excitations corresponds to multiply- 
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ing by a z on the red and a x r x on the green bonds. These 
combinations can be readily shown to commute with the 
Hamiltonian along the length of the string, the only non- 
commutatitivy being associated with the ends. Hence 
the energy cost to create the point defects at the surface 
is finite even for large separation. 



incurred during this process encodes the mutual statis- 
tics of this pair of particles. The e particles in Fig. [8] 
are created by Sf 3 = g\g 3 g\t% while the m particles 



are created by S2 

cmce Grace 

^24^13^24^13 — 



• = Tfrfrf . It is readily seen that 
/, demonstrating that e and m have 
mutual semionic statistics. 




FIG. 7. Deconfined fermionic excitations on the boundary. 
The path C12 shown in red connects a pair of vertices Vi and 
V2, each of which is associated with a fermionic excitation. 
The displaced path used to determine *Ci2 is indicated with 
a dashed blue line; edges that cross under this path are shown 
in green. Unlike in the bulk (Fig. [6]), on the surface the 
dashed blue line does not thread any plaquettes, and there is 
no string of defective plaquettes connecting the vertices. 

Statistics of Deconfined Surface Excitations: 



4 



FIG. 8. Mutual Statistics: The operation of braiding a 
pair of anyons (say the e and m particles) is captured by 
first creating a pair of e particles (red) followed by a pair of 
m particles (blue) in the manner shown. We now annihilate 
first the e and then the m particles to return to the vacuum, 
and examine the resulting phase. 



the boundary) have exactly the mutual statistics of the 
fermions e, m, and e respectively in our three fermion 
model. The crucial point is that if C12 and C34 are paths 
that cross once, the operators S^ 34 and S£ i2 (/i ^ v) do 
not commute. Consider, for example, the mutual statis- 
tics of the e and m particle. This can be obtained by 
first creating a pair of e particles, and then a pair of m 
particles as shown in Figure [8j and next destroying first 
the e particles and then the m particles. The net phase 




We now show that on the surface, the excitations cre- 
ated at the end-points of the string operators , S™ , S 

and Sq in Eq. (15) (where dj is an open curve in S 



FIG. 9. Fermionic statistics: The sequence of operations 
described in the text to exchange a pair of fermions. The 
positions of the fermions are indicated by the red dots. (The 
open red circles in (b) - (d) indicate the original positions of 
the fermions). Solid red arrows indicate the edge acted on 
by a string operator to move the fermion at this step; the 
solid red lines show where the string operators have acted at 
previous steps. The edges crossed by the dashed blue line are 
in *C: there is a phase of —1 every time a dashed blue line 
crosses a solid red line. 

Now consider the self statistics. The exchange of 
two anyons of the same type can be realized as shown 
in Fig{9] We begin with two anyons (a and b) of the 
same type at vertices i and i + x respectively. The first 
step in the exchange is to move the anyon a from i to 
i — y; then move anyon b from i + x to i — x; next 
move anyon a from i — y to i + x\ and finally move 
the anyon b from i — x to i. This process exchanges 
the two anyons. The whole procedure is realized by 
string operator .Sq i+£ ^c i+ , 

explicitly check that this operator is equal to 
fi = e,ra,£. For exmaple, when /1 = e. 



. We can 
-I for 



5 cv 



>b u 2> '3 ^ 



erg erf rf . 
-I. 



The total exchange string operator is 
Similar checks can be performed for 



then equal to 
\i = rn and e. 

Conclusion In this paper, we presented an exactly 
solvable model for the 3D bosonic SPT phase with 
time reversal symmetry (T 2 = 1). Our construction is 
based on the Walker- Wang prescription, which results 
in a topologically trivial 3D bulk state but a topologi- 
cally nontrivial 2D surface state that has three species 
of fermionic excitations with mutual semionic statistics. 
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Such a state cannot be realized in 2D without breaking 
time reversal symmetry, while our 3D model is explicitly 
time reversal symmetric. Therefore, the surface realizes 
a state that cannot be time-reversal invariant in a purely 
2D system, clearly indicating that the model has nontriv- 
ial SPT order in the bulk. The underlying principle of 
realizing a topologically ordered surface and a SRE bulk 
with a 3D loop gas can be generalized to other TQFTs, 
providing a natural platform for studying topologically 
ordered states arising at the surface of SPT phases. For 
example, the Z2 topological order can be naturally real- 
ized (see Appendix A). With appropriate transformation 
laws for excitations, this is predicted to be a possible 
surface termination for several other bosonic SPT phases 
with time reversal symmetry [17]. 
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APPENDIX 1: SURFACE Z 2 GAUGE THEORY 

Similar construction of 3D loop gas state can give rise 
to a SRE bulk state with 2D Z2 gauge theory on the 
surface. The only difference from the definition of the 3 
fermion model is the B matrix which is equal to the R 
matrix for the braiding of the electric quasiparticle e, the 
magnetic quasiparticle m and their bound state e 



u z 2 




(16) 



and eg = {el) T . 

The Hilbert space on each edge is again four dimen- 
sional with basis states labeled |l),|e),|m), and \e). If 
Pauli operators af (rf ) measure the Ising charges 
(Zr, 77 ^) flowing through the link z, then the exactly solv- 
able bulk Hamiltonian can be written as 



H 



J2a v -J2b p 



where 



with 



Av = -ai i o?+ii i 'Tf) 

b p = E M=e ,m,e HiedP S i 



(18) 




Open string operators are then given by 



s &» = n % n 



n 



(19) 



iec 12 je*c^er 



It can then be explicitly checked that the end of strings 
are confined in the bulk but deconfined on the surface 
and the surface anyons have the same statistics as the 
Z2 gauge theory in 2D. Of course, this theory can also 
appear in 2D with time reversal symmetry as defined. 
However, if anyons transform in more nontrivial ways 
this can lead to additional SPT phases [17]. This is left 
for future work. 
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